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Abstract. For the Euler equations governing compressible isentropic fluid flow with a 
barotropic equation of state (where pressure is a function only of the density) , local conser- 
vation laws in n > 1 spatial dimensions are fully classified in two primary cases of physical 
and analytical interest: (1) kinematic conserved densities that depend only on the fluid 
density and velocity, in addition to the time and space coordinates; (2) vorticity conserved 
densities that have an essential dependence on the curl of the fluid velocity. A main result of 
the classification in the kinematic case is that the only equation of state found to be distin- 
guished by admitting extra n-dimensional conserved integrals, apart from mass, momentum, 
energy, angular momentum and Galilean momentum (which are admitted for all equations 
of state), is the well-known polytropic equation of state with dimension-dependent exponent 
7 = 1 + 2/n. In the vorticity case, no distinguished equations of state are found to arise, and 
here the main result of the classification is that, in all even dimensions n > 2, a generalized 
version of Kelvin's two-dimensional circulation theorem is obtained for a general equation 
of state. 



1. Introduction 

Conservation laws and Hamiltonian structures are central to the mathematical study of 
fluid flow and have long been known for both the incompressible (ideal fluid) and compressible 
(inviscid fluid) Euler equations governing fluid flow in two and three dimensions. Over the 
past few decades there has been considerable mathematical interest in studying the Eulerian 
fluid equations in n dimensions [8]. 

One strong motivation came from the work of Arnold [6, 7J showing that the Euler equa- 
tions for incompressible fluids in n-dimensional spatial domains have an elegant geometric 
formulation as the geodesic equation on the Lie group of volume-preserving diffeomorphisms 
of the given domain of the fluid flow. This formulation gives an interesting geometrical sig- 
nificance to fluid conservation laws by interpreting them as geodesic first integrals related 
to invariance properties of the geodesic Lagrangian. Subsequently, the main group-theoretic 
aspects of Arnold's work were extended to the compressible Euler equations, first [131 HH] 
for isentropic fluids (whose entropy is constant throughout the fluid domain) in which the 
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pressure is specified to be a function only of density as given by an equation of state, then 
later [201 E] for adiabatic non-isentropic fluids (in which the entropy is conserved only along 
streamlines) where the pressure is given by a dynamical equation. 

The aim of the present paper and a sequel will be to give a complete picture of the 
conservation laws of kinematic type and vorticity type for general compressible fluids in n > 1 
dimensions for both isentropic and non-isentropic cases. By a kinematic conservation law we 
will mean a local continuity equation where the conserved density and flux depend only on the 
fluid velocity, pressure and density (but not their spatial derivatives), in addition to the time 
and space coordinates. Such conservation laws encompass the familiar physical continuity 
equations in two and three dimensions for mass, momentum and energy pTt [9l [12]. In 
contrast, a vorticity conservation law will refer to a local continuity equation for a conserved 
density and flux that have an essential dependence on the curl of the fluid velocity in a form 
exhibiting odd parity under spatial reflections. Examples of conservation laws with this form 
are helicity in three dimensions as well as circulation and enstrophy in two dimensions, which 
are well-known for incompressible fluid flow [TS], • 

To-date all of the known n-dimensional fluid flow conservation laws ^31 [13 HBJ belong to 
these two classes but have been derived through special methods that fall short of providing a 
complete classification. An interesting open question we will settle in this paper for isentropic 
compressible fluid flow is to find all particular equations of state for which the n-dimensional 
Eulerian fluid equations admit vorticity conservation laws or extra kinematic conservation 
laws. 

In section 2, as preliminaries, the general formulation of local continuity equations and 
integral conservation laws for the Euler equations for compressible isentropic fluids in n > 1 
dimensions is reviewed. In particular, we introduce necessary and sufficient determining 
equations for directly finding conserved densities of any specified form. By solving the 
determining equations for kinematic conserved densities, we obtain a complete classification 
showing that apart from mass, momentum, and energy, the only additional conservation laws 
of kinematic form consist of Galilean momentum (connected with center-of-mass motion) and 
angular momentum holding for any equation of state, plus dilational-type energies arising for 
polytropic equations of state where the pressure is proportional to a particular dimension- 
dependent power of the density. 

Next in section 3 we solve the determining equations to find all vorticity conservation laws, 
starting from the transport equation for the curl of the fiuid velocity. This classification 
yields an odd-dimensional generalization of helicity and an even-dimensional generalization 
of circulation and enstrophy, which are found to hold for any equation of state. We show 
that the generalized circulation has an equivalent formulation as a constant of the fiuid 
motion defined on the boundary of any spatial domain that is transported in the fluid. This 
new result gives a generalization of Kelvin's two-dimensional circulation theorem to all even 
dimensions n > 2. 

Finally, in section 4, we use the well-known Hamiltonian structure [23] of the compressible 
Euler equations to classify all Hamiltonian symmetries corresponding to the kinematic and 
vorticity conservation laws. In section 5 we make some concluding remarks, including a 
summary of our classification results stated in index notation. 
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A corresponding treatment of kinematic and vorticity conservation laws for adiabatic non- 
isentropic compressible fluids in n > 1 spatial dimensions will be given in a separate paper 

2. Compressible Isentropic Flows 

The Euler equations for compressible isentropic fluids in M" (in the absence of external 
forces) with velocity u(t, x) and density p{t, x) consist of 

9tu + u- Vu + ivp = 0, (2.1) 
P 

dtp+V-ipu) = 0, (2.2) 

with a barotropic equation of state for pressure 

p = P{p). (2.3) 

Throughout, we will use bold notation to denote vector or tensor variables and operators. A 
dot will denote the Euclidean inner product as well as stand for contraction between vectors 
and tensors, while a wedge will denote the antisymmetric outer product of vectors and/or 
antisymmetric tensors. 

Fluid conservation laws are described by a local continuity equation [211 SI IH] 

AT + ■ X = (2.4) 

holding for all formal solutions of f l2.ip -( l273l) . where T and X are some functions of t, x, p, u 
and their x-derivatives. Here Dt and A denote total time and space derivatives respectively. 
Physically, T is a conserved density with X being a corresponding spatial flux. In integral 
form, the continuity equation (12. 4p is equivalently described by 

^ / Td^'x = - [ X ■ M"~V (2.5) 
dt Jv Jdv 

where V is any spatial domain in M"' through which the fluid is flowing and n is the outward 
unit normal on the domain boundary dV. A physically more useful form for expressing 
fluid conservation laws 02.41) and (12.51) is obtained by considering a spatial domain V{t) that 
moves with the fluid [15]. Then the flux through the moving boundary dV{t) is ^ = X — Tu 
which is related to the conserved density T by the transport equation 

AT + u • AT = -(V • u)T - A ■ I (2.6) 

where D^ + u- A represents the total convective (material) derivative and V ■ u represents the 
expansion or contraction of an infinitesimal volume moving with the fluid. The corresponding 
integral form of a fluid conservation law in a moving domain is accordingly expressed as 

^ / Td"x = - / ^ ■ hd^-^a (2.7) 
dt Jv(t) Jdv{t) 

whereby jy^^^-^ Td^'x will be a constant of the fluid motion in V{t) C M"" if the net flux across 
dV{t) vanishes. 

The determining equations for finding conserved densities T are given by 



EpiVtT) = = E^iVtT) 
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(2.8) 



where Ep and are spatial Euler operators [1] with respect to p and u, and T>t is the total 
time derivative evaluated on solutions of the Euler equations fl2.ip - fl2.3p . (The explicit form 
of these operators is shown using index notation in section [2^21 ) These equations (I2.8P arise 
from the fact that spatial divergences Dy^ ■ X = —T>tT have a characterization pTl [TT] as 
functions of t, x, p, u and x-derivatives of p,u that are annihilated by both of the spatial 
Euler operators. 

A conservation law is locally trivial if the conserved density and spatial flux have the form 

T = D^-& = div@, X = ■ * - Dt@, (2.9) 

whereby the continuity equations fl2.4p and fl2.5p hold as identities for some vector function 
and antisymmetric tensor function ^ of t, x, p, u, and x-derivatives of p,u. The corresponding 
identity holding in a moving domain V{t) takes the form 

^ [ e nrf"-V = - / C ■ nrf"-V = - [ div^iTx (2.10) 
dt JdV{t) Jdvit) Jv{t) 

where ^ = —Dt@ — (Dx ■ 0)u is the spatial flux through the moving boundary dV(t). As 
it stands, fl2.10p has no physical content. However, an interesting observation is that if the 
moving-fiux ^ is divergence-free for all formal solutions of the Euler equations (I2.ip - (]2.3p 
then the quantity /gyj-j-, © ■ n(i"~^cr will be a nontrivial constant of motion. Namely, when 
u(i(:,x) and p(t,x) satisfy f l2.1l) -( l2l3i) . vector functions that satisfy the condition 

= -div$, = DtdivQ + ■ {{dw@)u) = ( A + u • + V ■ u)div@ (2.11) 

lead to nontrivial conservation laws of the form 

^ / 0-nrf"-V = O (2.12) 

d't JdV{t) 

for any boundary hypersurface dV{t) that is transported in the fluid. In particular, it is 
sufficient for to satisfy the condition of vanishing flux — ^ = Dt@ + [D^ ■ 0)u = for all 
formal solutions of (I2.ip - fl2.3p . We will call (12.121) a moving-boundary conservation law. 



2.1. Classification of kinematic conservation laws. We first consider kinematic con- 
servation laws as defined by the form 

T(t,x,p,u) (2.13) 

for the conserved density. In the case of polytropic equations of state, 

p = P{p) = Kp'^, K, J = const. (2-14) 

where pressure is proportional to a power of the density, all of the known local conserva- 
tion laws (12.131) of the compressible polytropic Euler equations fl2.1l) . fl2.2p . fl2.14l) in n > 1 
dimensions are summarized in the following table [1] 
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Conserved density T 


Description 


Number 


Exponent 


P 

pu 
pu A X 

p(tu — x) 

+ :pTP" = E 
p(| up + Klnp) 
tE — ipu ■ X 
t^E - tpu ■ X + ^p|xp 


Mass 
Momentum 
Angular momentum 
Galilean momentum 
Energy 
Energy 
Similarity energy 
Dilational energy 


1 

n 

n{n - l)/2 
n 
1 
1 
1 
1 


7 arbitrary 

)5 
)) 
)) 

7 = 1 
7 = 1 + 2/n 

)5 



We begin by stating a general classification of kinematic conservation laws with respect 
to all equations of state fl2.3p . Note that any conservation law of this form (12.131) is locally 
nontrivial since it does not contain x-derivatives of p or u. 

Theorem 2.1: (i) For a general equation of state fl2.3p . the fluid conservation laws (12.131) 
in any dimension n > 1 comprise a linear combination of mass, momentum, angular mo- 
mentum, Galilean momentum, and energy. In particular, for any spatial domain V{t) C 
transported in the fluid: 



_ / pdTx = 0, 

Vit) 

/ pu(i"x = 



d 
It 
d 
dt 
d 
di 
d 
di 
d 
di 



V{t) 

f 

V(t)' 
V(t) 
V(t) 



phd 



dV{t) 



pu A xci^a; 



px A hd"- V, 



dV{t) 



p{tu — x)d"'x 



tphd 



CT, 



dV(t) 

/ (ip|u|^ + pe)(i"x = - / pu-nrf""V 

Jv(t) JdV(t) 



where 



p = P(p), e= / p-'P{p)dp. 



(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 

(2.20) 



(ii) Modulo a constant shift in p, the only equation of state for which extra conservation 
laws (12.131) arise is the polytropic case (12.141) with dimension-dependent exponent 7 = 1 + -, 



p = np 



K = const. 



(2.21) 



The admitted conservation laws consist of a linear combination of a similarity energy and a 
dilational energy. In particular: 



^ / {tE- ipu ■ x)rf"x = - / p{tu - ix) ■ nd"-V, 
dt Jv{t) Jdv{t) 

f (t^E - tpu ■ X + |p|xp)(i"x = - /" pt{tu-x)-nd 

dt Jv(t) Jdv(t) 



CT, 



where 



E = lp\u\ + ^rihip 



\p\u\ + ^np 



(2.22) 
(2.23) 

(2.24) 



is the polytropic energy density. 

The proof of this classification theorem is given in section [22] using index notation. A sum- 
mary of the conservation laws fl2.15p - fl2.24p written in explicit component form is presented 
in section 5. 

We remark that these conservation laws were first derived [Hj for the case of irrotational 
fluid flow with a polytropic equation of state fl2.14p . The Euler equations for such fluids in 
n > 1 dimensions turn out to have a Lagrangian formulation when a velocity potential is 
introduced (i.e., V A u = implies u = V$), which allows local continuity equations to be 
classified in terms of point symmetries by means of Noether's Theorem [111121]. In particular 
mass conservation arises from invariance of the Euler- Lagrange fluid equations under shifts 
in the velocity potential. Invariance under space translations, rotations, Galilean boosts, and 
time translation respectively yields conservation of momentum, angular momentum, Galilean 
momentum, and energy. For the special polytropic equation of state fl2.21l) the conserved 
similarity energy arises from a particular combination of scaling and dilation invariance 
that produces a variational symmetry, while the dilational energy corresponds to an extra 
symmetry ^22j that is admitted only for this equation of state. 



2.2. Classification proof. The proof of Theorem 2.1 is based on explicitly solving the 
determining equations (12. 8p by tensorial index methods. We introduce the following index 
notation: x x*, u u*, Vp and Vu u\i (using a subscript comma to denote 

partial derivatives), and Di, where i = 1,2,. . . ,n] indices will be freely raised and 

lowered via the Kronecker symbols 6ij and 5*-' (which are components of the Euclidean metric 
tensor and its inverse on in Cartesian coordinates). The summation convention will apply 
to repeated indices. 
In index notation, 

ui = -v^v^, - p-'P\p)p:, pt = -{pu'),i, (2.25) 

are the Euler equations fl2.ip - fl2.3l) . The spatial Euler operators with respect to p and are 
given by 

For a conserved density of the kinematic form T{t,x''', p,u^), we have 

VtT = ~Tp{u'p^i + pu\i) - T^^iu^u'j + p-^P\p)p:) + Tt. (2.27) 
First applying the Euler operator Ep to (12.271) we get 

{kjpTpp - p-^P'{p)T^.^,)u\^ + p-^P'{p)T^^^^ + Tp,.u' + Ttp 

which is a linear inhomogeneous scalar expression in m*/ . Its coefficient must vanish, yielding 
the two equations 

6;,pTpp - p-^P'{p)T^^^, = 0, (2.28) 
p-'P\p)T^.,^ + Tp^.u' + Ttp = 0. (2.29) 

Next we apply the other Euler operator E^i to (12.271) . obtaining 

(5ijpTpp — p P' [p)Tuiu3) P ,^ — {pTuip — Tui)u\j + [pTuip — Tuj)u\i + pTpr^i -\- TuixjU-' + Ttui 



which is a hnear inhomogeneous expression in p -^, u\j and u\i. First we see the coefficient 
of yields the same terms as in fl2.28p . Next, since u\j and are hnearly independent 
in n > 1 dimensions, their coefficients must separately vanish, which yields 

pT^.p - = 0. (2.30) 

This leaves the inhomogeneous terms 

pTpxi + T^i^jvP + Tt^i = 0. (2.31) 

Hence the determining equations consist of fl2.28l) - fl2.31l) to be solved for T. We start from 
equation (12.301) . which is ffist-order linear in p. By integrating with respect to p, and then 
doing a trivial integration with respect to u*, we obtain 

T = pf{t,x\u')+g{t,x\p). (2.32) 

Substituting fl232D into (12:281) gives 

Sijpgpp = P'ip)fu'ui (2.33) 
which separates with respect to p, into two equations 

= c{t,x')5kj, gpp = c{t,x')p-'P'{p), (2.34) 
where c(t, x*) is a constant of separation. Integration of (I2.34p yields 

/ = ci(t, x') + Cj{t, x')u^ + |c(t, x')u^Uj, (2.35) 

g = Co{t, x') + C2{t, x')p + c(t, x')pe, (2.36) 

e = j p-'P{p)dp, (2.37) 

where Cq, Ci, C2 and Cj are constants of integration with respect to p and m*. Thus, we have 

T = Co + cp + Cipu^ + ciJ^pu^Ui + pe) (2.38) 

where c = ci + C2. Since the term cq is trivially conserved, we can put cq = 0. 

By substituting (I2.38P into fl2.29p and (12.310 . in each case we get a cubic polynomial in 
terms of whose separate coefficients must vanish. This leads to the following system of 
equations: 

c^ = 0, (2.39) 

Cj,i + Cij + Ct5,j = 0, (2.40) 

du + c , = 0, (2.41) 

ct + P'(p)q; + (pe)'Q = 0. (2.42) 

To proceed, we note (]2.39p immediately implies c = c{t). Then (I2.40p has the form of a 
time-dependent dilational Killing vector equation on q. To derive the solution, we ffist take 
the antisymmetrized derivative of (12.401) . i.e. differentiating with respect to x'' followed by 
antisymmetrizing in j and k, which yields 

{5j,k - 4,i),* = 0. (2.43) 

Similarly, by taking the curl of (12.410 . i.e. differentiating with respect to x^ and antisym- 
metrizing in i and j, we obtain 

(Q,i - 5j,i)t = 0. (2.44) 



Hence fl2:i3|) and flCTj) give 

^hj ~ ~ '^^lij (2-45) 

where Cuj is an antisymmetric constant tensor. Adding (12.451) to fl2.40p . we get 

25,,, = 2Cuj - c'{t)5i, (2.46) 
and thus, by integration with respect to , 

Q = Co^{t) + CujX^ - lc{t)xi, Cuj = -Ciji. (2.47) 

Then (12.411) becomes 

c^ = -C'oi + \d'xi (2.48) 

which yields 

a = C2(t) - CoiX* + Yx'xi. (2.49) 
Finally, from (I2.42p . by using (12.491) and the trace of (I2.40p . we get 

C'2 - C'o^x' + lc"'x% + hlnP' + {pe)'y = 0. (2.50) 
Splitting (12.501) with respect to x^ yields 

C;', = 0, c'" = 0, (2.51) 

C2 + ((pe)' - ^nPy = 0. (2.52) 

From (I2.5ip we have 

Coi = CLoi + ttiit, c = Bq + bit + h2t^ (2.53) 

with constants aoi, Oi,, 60 > ^i; ^2- Differentiating (I2.52p with respect to p gives 

((pe)' - \nP')'d = (2.54) 

which leads to the following two cases. 

Case c' = 0: Hence c = bo and 61 = 62 = 0, which implies C2 = const, from (12.521) . Then 
(^M) and (nmf yield 

Ci = aoi + aut + Cujx\ c = C2 - aux\ (2.55) 
whence from (I2.38P we obtain 

T = C2P + aoipu" + aup{tu' - x*) + Cajpu^'x^ + bopHu'ui + e) 

where C2, Cuj = —Ciji, a^i, au, 60 are arbitrary constants and e is a function of p given by 
(E3ZD. 

Case (pe)" - f P" = 0: Simplifying (pe)" = pe" + 2e' = p-^P' through (^M), we get a 
linear ODE 

npP" - 2P' = (2.56) 

which has the general solution P(p) = do + rfip^"^^/", where do, di are constants. Note that 
we can put do = since the pressure p = P{p) can be shifted (without loss of generality) by 
an arbitrary constant. This implies 

P = dipi+2/", e = indip^/" + d2 (2.57) 

from (I2.37p . with a constant of integration d2. Thus, (I2.52p becomes + d2c' = whence 

C2 = 63 - d2c{t) = 63 - ^260 - d2bit - d2b2t^ (2.58) 
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where 63 is an integration constant. Then (12.470 and (12.490 yield 

Ci = aoi + aiit + CiijX^ — + 2b2t)xi, c = C2(t) — aux^ + jb2X^Xi. (2.59) 

Substituting (12.581) and (12.591) into (I2.38p . we find that the terms involving d2c(t) cancel out, 
giving 

T = 63P + aQipu^ + aiipitu^ — a;*) + Cujpu^x^ + bop{^u^Ui + e) 
+bip{lu'{tui - Xi) + te) + b2p{\{tu' - x'){tui - Xi) + t^e) 

with 

e = InP/p, P = rfip^+"/2 (2.60) 

where di, bo, bi, 62, &3, aoi, clu, Cuj = "Ciji are arbitrary constants. 
This completes the proof of Theorem 2.1. 

3. VORTICITY CONSERVATION LAWS 

In n > 1 dimensions, the curl of the fiuid velocity is the antisymmetric tensor 

a; = VAu (3.1) 

which satisfies the identities V A uj = and V • = Au — V(V ■ u). There is a natural 
odd-parity expression that can be constructed purely out of products of u) and the spatial 
orientation tensor e as follows. (Recall, e is a rank n totally skew-symmetric tensor whose 
components in Cartesian coordinates for R" are given by the Levi-Civita symbol. In partic- 
ular, up to a choice of sign, e is determined by its properties Ve = and |ep = e- e = n!.) 
When the spatial dimension is even, say n = 2m, 

ti7 = e-((VAu) A---A(VAu)) (3.2) 

^ v ' 

m times 

defines a vorticity scalar that has odd parity since e changes sign under spatial refiections. 
Similarly, when the spatial dimension is odd, say n = 2m + 1, the analogous expression 

z:c7 = e-((VAu) A---A(VAu)) (3.3) 

^ V ' 

m times 

defines a vorticity vector with odd parity under spatial refiections. These expressions (13.21) 
and (13. 3p can be written in the more compact notation *{uj'^), where * is the Hodge dual 
operator (acting by contraction with respect to e) and m = [n/2] is a positive integer. 
We now define vorticity conservation laws to have the form 

T(p,u,*(a;-)) (3.4) 

for the conserved density, such that the expression (13. 4p possesses odd parity under spatial 
refiections as follows. Let be a refiection operator defined with respect to a spatial unit 
vector £ in R". Specifically, reverses all vectors parallel to £ while leaving invariant all 
vectors in the hyperplane orthogonal to £ through the origin. Note extends to act on 
tensors by multi-linearity and acts as the identity on scalars. Then for any choice of £, Vi 
satisfies the properties V^e = —e and Vj = I, so consequently the parity of a conserved 
density (13. 4p will be odd under spatial refiections iff 

T{Vp, Vu, *{Vuj^)) = -VTip, u, *(a;™)). (3.5) 
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Note that, in contrast, the parity of all the kinematic conserved densities fl2.15l) - fl2.23l) is 
even, 

T(t, Px, Vp, VvL) = VT{t, X, p, u). (3.6) 

Our main result will now be a complete classification of vorticity conservation laws for all 
equations of state fl2.3p . 

Theorem 3.1: For any equation of state fl2.3p the only nontrivial fiuid conservation laws 
f l3.4p - fl3.5l) in dimensions n > 1 are given by helicity 

u-vjrx = - [ (/(p) - i|un-cc7-nrf"^V (3.7) 

dt Jv{t) Jdv{t) 

for odd dimensions n = 2m + 1, where 

f{p) = [ p-'P'{p)dp = - + e, (3.8) 
J P 

and generalized enstrophy 

^ / pf{w/p)rx = (3.9) 
"'^ Jv{t) 

for even dimensions n = 2m, where / is any nonlinear odd function of w/p. In particular, 
there are no special equations of state that admit extra vorticity conservation laws. 

The helicity and enstrophy conservation laws were first derived by means of a Hamiltonian 
Casimir analysis fiOi [8] , as we will discuss further in section HJ which is a more restrictive 
analysis than directly solving the determining equations fl2.8p for conserved densities. By 
comparison, our classification has more generality and actually holds without imposing the 
odd-parity condition fl3.5p if we consider conserved densities 03.40 that just have an essential 
dependence on uj. The proof is given in section 3.1 and a summary of the conservation laws 
fl3.7p - fl3.9l) in explicit component form is shown in section 5. 

The local continuity equations underlying these conservation laws fl3.7l) - fl3.9l) are readily 
obtained from the transport equations satisfied by zcr and vu. Through (13.11) . we first note 
the identity 

u ■ = u • Vu — |V(u ■ u) 
whence the Euler equation (12.11) for the fiuid velocity can be written in the form 

9tu + u-a; + V(||up + /(p)) = (3.10) 

where f{p) is given by (13. 8p through the equation of state (12.31) . Taking the curl of (13.101) 

we get 

dtUJ = V Aiuj -u). (3.11) 

Substitution of (13.111) into the time derivatives of (13.21) and (13.31) then yields the respective 
transport equations 

dtw + u • Vts7 = —(V • u)ti7, dt^dj + u ■ V-ccT = Tu • Vu — (V ■ u)^^?, (3.12) 
or equivalently, 

dtvj = -V ■ (wn), (3.13) 
dt-ca = V ■ {vj Au). (3.14) 
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The vorticity transport equations (13.131) and (13.141) each have the form of a local continuity 
equation holding in even (n = 2m) and odd (n = 2m + 1) dimensions, with respective 
conserved densities T = w and T = vj. As we will now show, both of these vorticity 
conservation laws are locally trivial and therefore fall outside of our classification theorem. 
From the explicit expressions (13. 2p and (13. 3p for the vorticity densities in terms of the curl 
of the fluid velocity lj = V A u, we find that in both even and odd dimensions 

T = >K(a;™) = V- *(u Aa;™-^) = rfii;0 (3.15) 

is a spatial divergence. The fluid equations (I3.10p and (13.111) yield the time derivative of 
e = *(u Aa;'"-!) to be 

dt@ = -viA*{uj'"') - W --^ (3.16) 

where ^ = (||up — /(p)) * (cj™^^) + u A (u ■ *(cl>'"^^)) is an antisymmetric tensor. Since in 
both even and odd dimensions the expression 

X = uA*(a;™) (3.17) 

is the spatial flux, X = wu and X = u A tu, arising from the vorticity transport equations 
(I3.13P and (13.140 . we then see (I3.17P and (13.150 have the form (12.90 of a locally trivial spatial 
flux and a locally trivial conserved density. 

It is interesting to investigate the corresponding vorticity flux ^ = X — Tu through a 
moving domain boundary in the fluid. 

In odd dimensions n = 2m + 1, we can write 

vu = -V ■ {u-W) (3.18) 

where 

W = *(a;'"-i) = ^9^117 (3.19) 

is a totally skew-symmetric tensor of rank 3. Hence for any domain V{t) transported in the 
fluid, we obtain 

- ^ / vjd^'^+^x = ^ [ W ■ (u A n)rf2™(T = [ u{vj ■ h)d^"'a (3.20) 
dt Jv{t) dt Jgv{t) JdV{t) 

for all formal solutions of the dynamical equations (I3.10p and (I3.14p . Since the moving- flux 
through dV{t) fails to vanish, (13.200 does not yield a constant of motion. 
For even dimensions n = 2m, we have 

zu = -V ■ (u-w) = V ■ {w ■ u) (3.21) 

where 

w = *(a;'"-i) = i9^tJ7 (3.22) 

is an antisymmetric tensor. Then for all formal solutions of the dynamical equations (I3.1O0 
and (I3.13p . we obtain 

^ zu(f"'x = ^ [ w (u A n)d'™- V = 0. (3.23) 



dt Jyi^t) dt jQV(t) 

Thus (I3.23P yields a constant of motion for the moving boundary dV{t) in M^™. It describes 
an even-dimensional generalization of Kelvin's circulation for isentropic fluid flow in two 
dimensions [9j. Specifically, when n = 2 (m = 1), the vorticity scalar is given by 

uj = e- uj = V ■ (*u) (3.24) 
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and hence w = e is the spatial orientation tensor. This yields w ■ (u A n) = u ■ *n, and then 
the conservation law fl3.23p becomes Helmholtz's circulation theorem [12\ 

d f d f 

zu(fx = -J- u ■ *nda = (3.25) 



dt Jy(t) dt jQVit) 

where dV{t) is a closed curve in M^, *n is a unit tangent vector along dV{t) and da is the 
arclength element. By writing ds = *hda we see fl3.25p states that the line integral 

u ■ ds (3.26) 

defining the circulation of the fiuid velocity around a curve transported in the fiuid is a 
constant of the fiuid motion. For higher dimensions n = 2m {m > 1), we can write the 
moving-boundary conservation law fl3.23p in an analogous form by noting 

w- (uAn) = (uAo;™-^) ■ *n (3.27) 

where *n is the volume tensor for the boundary hypersurface dV{t) in M^™. Then (13.231) 
becomes 

^ md^'^x = ^ [ (u A cj'"-^) ■ *nrf2"^-V = (3.28) 



which states that the (hyper) surface integral 

(uAo;"^-^) -(iA (3.29) 



'dV{t) 

is a constant of the fiuid motion, where dA = *hd'^"^~^a denotes the volume element for the 
moving-boundary (hyper) surface dV(t) in even dimensions n = 2m analogous to ds = *hda 
for moving-boundary curves in two dimensions. 

Proposition 3.2: The only moving-boundary conservation law (12.121) of vorticity type is 
the generalized circulation (13.280 holding for any equation of state in all even dimensions. 

The proof of this classification is given in the next section. 

3.1. Classification proof. We will use the same index notation introduced in section 2.2 
for the proof of Theorem 2.1. To begin, we write out the component form of the fiuid curl, 
the vorticity scalar and vector, along with their transport equations: 

w^w = e'^^'-'^^^Ui,,, ■ ■ -o^i^j-^, m = n/2, 
^^w' = e'^'^'-'-^''-^u,,k, ■ ■ ■ m = (n - l)/2. 



Wt = —{wu^)^i = —ZUu\i — U^TUi, 

In addition, we will need the component form for 

w ^ w'^ = e'^''^'-'"^-'^^~'uJi^j, ■ ■ ■ uji^_,j^_„ m = n/2, 

W ^ W'^' = e'^''^^''^-^-~^'-^u;,,k, ■ ■ ■ m = in- l)/2, 
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Here e*^'"*" = el'i'"*"! ^ e, Sij = 5(ij) ^ g are the components of the spatial orientation tensor 
and the Euchdean metric tensor; round brackets denote symmetrization of the enclosed 
indices, and square brackets denote antisymmetrization. 

The proof of Theorem 3.1 and Proposition 3.2 proceeds by explicitly solving the deter- 
mining equations (12.81) for conserved densities of vorticity type in even and odd dimensions 
n > 1. Recall, the Euler equations fl2.ip - fl2.3l) are given by 

ui = -u^u\,-f'{p)p:, pt = -{pu')^„ (3.30) 

where 

f{p) = P\p)/p ^ 0. (3.31) 
Case n = 2m: For a conserved density of the form T(p, w), its time derivative is given 

by 

-VtT = {pu'),T^ + {u'u\,+p:f')T,. + {wu^),T^ 

= u\p,Tp + u\iT^, + ro^T^) + u\i{pTp + wT^) + p^f'T^.. (3.32) 

In (I3.32p . we see the coefficient of equals DiT by the chain rule, which allows us to write 
the corresponding terms as u^DiT = Di{u^T) — u\iT. Hence 

AT = D,{-u'T) + u\,A - p^f'T^. (3.33) 

where 

A = T-pTp- wT^. (3.34) 
To begin we substitute f l3.33p into the first determining equation 

= EpiVtT) = tiuAp + v'^T^^ + w/T,,^ (3.35) 
where we have introduced the notation 

iiu = 6iju'^ = u\i, u'^ = i(M*/ + u{') = u^\^\ (3.36) 
Since T does not contain any derivatives of cu, the coefficient of in (13.350 must vanish, 

Tu^^ = 0. (3.37) 
Integrating fl3.37p and dropping a kinematic term c(p, u") that does not involve we get 

T = a{p,to). (3.38) 
Then (13.350 reduces to = tiuAp and hence we obtain 

= Ap = -papp - wa^p. (3.39) 
This is a first-order linear PDE for a^, which yields 

a = h{w) + pc{w/ p). (3.40) 
Thus from ( 13.380 and ( 13.400 we have 

T = h{w) + pc{w/p) (3.41) 

and 

VtT = Di{-u'T) + u\,{h-vjh'). (3.42) 
Now we substitute (I3.42p into the second determining equation 

= Eui{T>tT) = zuiB + mtTu^^zujiB + mvj^-'wjiiiuB' (3.43) 
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where 

B = zub". (3.44) 

Note that b{w) has no dependence on trw since w does not contain u^^ (which is hnearly 
independent of Wij). Consequently, the coefficients of all terms in fl3.43p must vanish, yielding 
B = B' = so thus b" = 0. Hence, we obtain 

b = bo + bizu (3.45) 

with constants 60, ^i- Therefore, (13.451) and (13.411) give the result 

T = biw + pc{w/p) + bo (3.46) 

with 

VtT = -Di{u'T) (3.47) 

from (I3.42p . where bi, bo are arbitrary constants and c is an arbitrary function oiw/p. Since 
the term 60 is trivially conserved, we can put 60 = 0. 

We now note that the ffist term in (13.461) is a trivial conserved density 

biw = bi{w'^Uj)^i = Di{biw'^Uj). (3.48) 

Its corresponding moving-flux vanishes, ^* = 0, due to the form of (I3.47p . Similarly, if 
c = cw/p is a linear function of w/p, where c is a constant, then the second term in 
(I3.46P becomes a trivial conserved density, pc{w/p) = cm. This proves Theorem 3.1 and 



Proposition 3.2 in the even-dimensional case. 

Case n = 2m + 1: Here we will need the identities 

w\i = 0, (3.50) 

w'uij = 0, (3.51) 

which are consequences of 

W'^'' = W^'^''\ (3.52) 

W'^'^uj^^^k = W'^^Uj,k = 0, (3.53) 

iOjk-^^^eij^kvjmkra = uluj j[iUJ j^kj " " " ^j^k,„] = n\uj[jiUJ j^kl ■ ■ ■ ^jrakm] = 0, (3.54) 



where (13.541) follows from the fact that there are no totally skew-symmetric tensors of rank 
2m + 2 > n = 2m + 1. 

Now for a conserved density T(p, m*, cu'), by the same steps followed in the previous case 
to evaluate the time derivative, we have 

VtT = Di{-u'T) + u\i{T - pTp - w^T^,) + u\jW^T^, - p^f'T^,. (3.55) 

To proceed we substitute (13.550 into the ffist determining equation = EpiVtT). This 
yields the terms 

= u'^WjTp^.. - iiu{pTpp + w'Tp^^) + u'^{wjTp^^ + f'T^^uo) + ^^'fTu^^j (3.56) 

using the notation (I3.36p . To start, we observe the ffist term in (I3.56P vanishes by the 
identity (13.510 . Next, since T does not contain any derivatives of cu*, the last term in (13.560 
must vanish modulo the identity (13.500 . This implies 

T^^^j = a{p,u'',w'')6ij. (3.57) 
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Applying the derivative operator to (13.571) and antisymmetrizing in [j/c], we get 

dija^k — SikttuJ = 0. (3.58) 

The trace of this equation with respect to {ij) yields {n — l)auk = 0, whence in n > 1 
dimensions we obtain 

a^k = 0. (3.59) 
By also applying the derivative operator d^k to (13.571) . we similarly get 

a^k = 0. (3.60) 

Integration of (IXFTD . (IX^ . (IXBUD then yields 

T = a{pywi + b{p,w'). (3.61) 

Here we have dropped an integration constant c(p, m*) since it does not involve ccr* (i.e. it is 
of kinematic form). Hence the determining equation (13.561) reduces to 

= u'^{Bij -6i,A) (3.62) 

with coefficients 

A = phpp + w%p, (3.63) 

Bij = zu(^jbi)p + a'w(^jUi)^ (3.64) 

in terms of the notation hi = b^i. Since tu* does not contain m*-' we see that b{p,zu^) has 
no dependence on m*-' (or tiu). Consequently, (13.621) implies that the coefficient of u^^ must 
vanish, 

Bij = A5ij. (3.65) 

By considering the product of (13.641) with Wk^i and antisymmetrizing in [il] and [jk], 
we find wiiBi][jtUk] = 0. The same antisymmetric product applied to (13.651) then implies 
Awi^i6i][jZUk] = which gives 

A = (3.66) 

and hence 

Bi^ = 0. (3.67) 
Now by taking the product of (I3.67P with zuk antisymmetrized in [jk], we get 

TU[kbj]p + a'zu[kUj] = 0. (3.68) 
This can hold only if a' = and bjp = zujc{p, w^). Then (I3.64p becomes 

Bij = czUiWj (3.69) 
whence (I3.67P implies c = and so bjp = 0. Thus we have 

b^^p = (3.70) 

whose solution is 6 = 6(ti7*) +c(p). Since c(p) contributes only a kinematic term in (I3.6ip . it 
will be dropped hereafter. Then 

b = 6(tJ7*), a = const. (3-71) 

satisfies both (I3.67P and (13.661) . 

Thus we have, from (13.711) and (13.611) . 

T = au'wi + b{w') (3.72) 
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and from fl3.55p . 

VtT = Di{-u'T) + u\j(h - ui%^) + u\jUJ^{aUi + h^^) - p,iaf'zu\ (3.73) 

Through (13.501) . we note p^iafw"^ = Di{fzu'') and u\jZu^Ui = DjiJ^u^UiW^) while u\jzu^b^i = 
{v}^ + u}^^)wjh^i = u^^Wjh^i by (I3.5ip . Thus we have 

VtT = Di{a{\u^Uj - f)w' - u'T) + u\j{h - UJ%^) + u'^zujb^^. (3.74) 

The second determining equation = E^i {T>tT) thereby yields the following terms 

= A,w\' + B'h^\j + mW'^'itrujAk + u<^\jB,hk + w\j{tiuAk^i + m^'^^.^.^^O) (3-75) 
where 

Ai = -w^bji, Bijk = Sk{ibj) + W(^ibj)k, (3.76) 
with the notation 6j = b^i, bij = b^i^j. Using the identities 

= W^pguP'^j, u^\' = u'^^ + fcu^'(^'\ trw/ = M^'^fc - luo^\k, 
and collecting like terms in (I3.75p . we get a linear homogeneous expression in u^^\ uj'^'^\^\ 

= u^'^'W'^^.iBkiH + SkiAh) - IJ^'W^'.^iBkih + SmAh) + uj^\m\i{B,,u + 5,.v4,) 

+mu'''uP\m%gW'',,{BuH^, + 5kiAf,^,). (3.77) 

Since uj''^''p is linearly independent of m-''^' in n > 1 dimensions, their coefficients in (I3.77P 
must vanish. From the coefficient of to''^'- ^^^u^'^, we have 

W'pgW',,iBgi,i^. + 6ghAi^.) = (3.78) 

which yields 

= 0, A^, = 0. (3.79) 

These two equations imply Bkji must be a constant tensor and Ai must be a constant vector. 
Next, the coefficients of u^''^ and uj^^\^'> yield the equations 

+ = o, (3.80) 

W\h{Bijk + 5ijAk) - 2W\i^j{Bi)gk + 5i)gAk) = 0. (3.81) 

We note that antisymmetrizing (13.800 in [ij] leads to (I3.8ip after indices are renamed. Hence 
only f l3.80p needs to be considered. Taking the trace of f l3.80p in (ij) and using the identity 
dSJni), we obtain 

W'^^kBDjH = 0. (3.82) 

Such an algebraic equation can hold only as a consequence of the skew-symmetry property 
fl3.52p . so fl3.82p is satisfied only if Bij^ = B(^ij)h is of the form 

Bijh = ^h{iCj) + SijCh + Cijh (3.83) 

where q, Ch are constant vectors and c^h = C(^ijh) is a trace-free totally symmetric constant 
tensor. Hence fl3.80p reduces to 

W\y5ki) (4 + Ah) + W\(^jCki)h = 0. (3.84) 
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The trace of fl3.84p in (kl) leads to Ch + Ah = (since Chk'' = and W^^iCjkh = 0) which 
imphes 

W\^jcu)h = 0. (3.85) 

By the same argument that led to (13.831) . the only totally symmetric tensor that can satisfy 
(I3.85P is Ckih = S(kiCh), which is trace-free only if Ch = 0. Thus, we have 

Bijh + SijAh = 5h(iCj), Cj = const. (3.86) 

giving the solution of equation (I3.77p . Substitution of expressions (13.761) into (13.861) then 
gives us 

5k(ibj) + zu(ibj)k = h(iCj) + dijw'-hik. (3.87) 
Antisymmetrizing (13.871) in [jk] yields 

8i[khj] = W[jhk]i (3.88) 

where hj = cj — bj — 2w^bkj. Taking the product of (13.880 with wi antisymmetrized in [jkl], 
followed by taking the trace in (ik), we get by'coij = 0. This implies 

bj = d{w')wj (3.89) 



whence (I3.88p yields 

bki = d{w^)5ki + d{w^)wk'Wi (3.90) 

and thus 

Cj - bj = {3d + 2dw^Wk)-ojj- (3.91) 
Substituting fl3:90l) and fl3:9B into (13:871) . we get 

dzuiZUjZUk = 3(d + dw^wi)5{ijWk) (3.92) 

which directly implies 

= J = 0. (3.93) 

From (I3.90p we thus have b^k^i = bki = and so, by direct integration, 

b = bk'OJ^ , bk = const. (3.94) 

where we have dropped a constant term since it does not involve . 
As a result, from (I3.72p we have 

T = au'wi + biw' (3.95) 

with arbitrary constants a, 6j, while from (I3.74p 

VT = -Di{a{f - \u^Uj) + u'T - bju'w') (3.96) 

where the last term comes from 

6j («%'),, = bjw'u\i = bjWiiu^' + uj^') = bjWiu'^ (3.97) 

due to (I330|) and (133TD . 

To conclude the proof of Theorem 3.1 and Proposition 3.2, we note that in (I3.95P the term 

b,w' = Uw'^^'uk),, = D^ihW^^'uk) (3.98) 
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is a trivial conserved density, while its corresponding moving-flux from fl3.96p is given by 
^* = —bjU^w^, which fails to be divergence-free, 

Di^' = -bj{u^w')^i = -%w'u\, = -hjWiu'^ ^ 0. (3.99) 

4. Correspondence between conserved densities and Hamiltonian 

symmetries 

The n-dimensional compressible Euler equations fl2.ip - fl2.3l) have the well-known Hamil- 
tonian formulation [231 



u 



5E/5u\ __/p-i(VAu)- -V 



where e = J p^^P{p)dp is the internal energy density and Ti. is called a Hamiltonian operator. 
This means Ti determines a Poisson bracket [21] 



{J'^GU = J{6F/6n 5F/5p)n(^fJ^^^^drx (4.2) 

= j p-^ ( V A u) ■ {5F/5VL A 5G/5u) + 5G/5vl ■ V5F/5p - 5F/5u ■ V5G/5p (Tx 

having the properties that (modulo divergence terms) it is antisymmetric and obeys the 
Jacobi identity, for arbitrary functionals T = j FdJ^x and Q = J GdP'x where F and G 
are functions of t, x, p, u and their x-derivatives. Here 5/5vi and 5/5p denote variational 
derivatives, which coincide with the spatial Euler operators -Eu and Ep when acting on 
functions that do not contain time derivatives of p and u. 

To check the Hamiltonian structure (14.11) produces fl2.ip - fl2.3l) . we note 6E/6u = pu and 
SE/Sp = ||up + e + p~^P{p), which yields 

p-^{V A u) ■ 6E/6u = (V A u) ■ u = |u| V|u| - u ■ Vu, 
-V ■6E/6u = -V ■ (pu), 
-V6E/6p = -iV|u|2 - V(e + p-^P(p)) = -|u|V|u| - p-^P'{p)Vp, 
whence we obtain 

6E/6u \ _ / p-i(V A u) ■ 6E/6u - V6E/6p 
SE/6p J ~ \ -V ■ 6E/6u 

-u- Vu-p-iVP(p) \ ^ f dtu 
-V ■ (pu) J [dtp 



n 



(4.3) 



Through this formulation the Hamiltonian operator Ti gives rise to an explicit mapping 
that produces symmetries of the compressible Euler equations from conservation laws as 
follows. Recall, fluid symmetries [15] are described by an infinitesimal transformation 

Xu = r), Xp = fj, (4.4) 

on all formal solutions of fl2.ip - fl2.3p . where t) and fj are some functions of t,x, p, u, and 
x-derivatives of p, u determined by infinitesimal invariance [2T1 [TU] of the Euler equations 

Dtf) + u ■ Vr) + r) ■ Vu + V {p^^ P' {p)fi) = 0, Dtf] + V ■ (r)u + pfj) = 0. 
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Now if T is a conserved density of the Euler equations fl2.1l) - fl2.3p then the mapping 

5T/5vi \ _ Y ( ^ 



can be shown (cf. general results in [21]) to yield a symmetry X = r)J9u + f^dp, given by 

r) = -p-i ( V A u) • 5T/5u + V6T/6p, r) = V ■ 6T/6u. (4.6) 

In particular, as seen from (I4.3p . the conserved energy density fl2.19p yields a time translation 
symmetry 

T = ip|up + pe^Xu = -dtu, Xp = -dtp. (4.7) 

For the other kinematic conserved densities listed in Theorem 2.1, we find the following 
correspondences: mass density fl2.15p yields a trivial symmetry 

T = p ^ Xu = = Xp; (4.8) 

momentum densities (12.160 yield space translation symmetries 

T = pu ^ Xu = Vu, Xp = Vp; (4.9) 

angular momentum densities (12.170 yield rotation symmetries 

T = puAx^Xu = (xAV)u-uAg, Xp = x A Vp; (4.10) 

Galilean momentum densities (12.180 give rise to Galilean boost symmetries 

T = p(tu-x) ^Xu = tVu-g, Xp = tVp. (4.11) 

Here g is Euclidean metric tensor on M" (recall, in Cartesian coordinates, the components 
of g are given by the Kronecker symbol). 

In the case of a polytropic equation of state (12.140 with special exponent 7 = 1 + 2/^, the 
similarity energy (I2.22p yields 

T = tE- |p(u ■ x) ^ Xu = -{tdtu + |(x ■ V)u + |u), 

Xp = -(t9tp + |(x-V)p+|np) (4.12) 

which is a scaling (similarity) symmetry; and the dilational energy (12.230 yields 

T = t^E - tp(u ■ x) + ip|xp ^ Xu = -{t^dtVL + t(x ■ V)u + tu - x), 

Xp = -{t^dtp + t{^-V)p + ntp) (4.13) 

which we call a Galilean dilation symmetry because it preserves x — tu and x/t. Here E is 
polytropic energy density (12.240 . 

In contrast, for the helicity and generahzed enstrophy densities (l3.70 - (l3.9p . we get: 

T = u--cc7^Xu = = Xp; (4.14) 

T = pf{w/p) ^ X\i = Q = Xp. (4.15) 

Thus we have the following classification result. 

Proposition 4.1: In all dimensions n > 1, the nontrivial infinitesimal symmetries pro- 
duced from the kinematic conserved densities (I2.15p - (l2.23p under the Hamiltonian mapping 
(14. 5 p consist of space translations (14. 9p . rotations (I4.1O0 . Galilean boosts (14.110 . and a time 
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translation (14. 7p for general equations of state, plus a similarity scaling (14.120 and a Galilean 
dilation (I4.13P for the special polytropic equation of state (12.211) . 

All of these Hamiltonian symmetries (I4.7I) - (I4.15I) can be seen to have the form of infini- 
tesimal point transformations 

f] = T] — rdtu — ^ ■ Vu, fj = T] — rdtp — ^ ■ Vp (4.16) 

so thus 

X = fi\du + fidp^X = Tdt + ^J9x + rildu + vdp (4.17) 

where rj, rj are functions of t, x, u, p, while r, ^ are functions only of t, x. These symmetries 
have the following geometrical description (proven in section 5). 

Proposition 4.2: Let C(x) be any solution of the dilational Killing vector equation 
£(^g = Qg, Q = const., on R", and let x(x) be any irrotational solution of the same equation, 
C-)^g = Qg, Q = const., V A X = 0, on M". Then the Hamiltonian symmetries (I4.8I) - (I4.13I) 
corresponding to the kinematic conserved densities (12 .161) . (pTT7|) . (12 . 1 8"|) . (12 . 22p . (12 . 231) have the 
form 

X = C\d^ + ntdt + \nnpdp + (^]u + \{V a C) ■ u)J9u, (4.18) 
X = tx\d^ + lnt% + \nntpdp + {\tnvL - x)\du- (4.19) 

A comparison of these Hamiltonian symmetries with all of the well-known point symme- 
tries [221 [IS] admitted by the compressible Euler equations (I2.ip - (l2.3p in dimensions n > 1 
gives the following classification results. 



For general equations of state (12.30 : 



Point Symmetry X 


Description 


Number 


Hamiltonian 
Gorrespondence 


X A 9x + u A 9u 
td^ + 9u 
tdt + xJ 9x 


Time translation 
Space translations 
Rotations 
Galilean boost 
Dilation 


1 

n 

n{n - l)/2 
n 
1 


Energy 
Momentum 
Angular momentum 
Galilean momentum 
Nil 



For polytropic equations of state (I2.14p : 



Point Symmetry X 


Description 


Exponent 


Hamiltonian 
Gorrespondence 


xJ 9x + uj 9u + Tpipdp 
tdt + ^xj^x - luj^u - fp9p 
t^dt + txj 9x + (x — tu)J d-a — ntpdp 


Scaling 
Similarity Scaling 
Galilean dilation 


7^1 
7 = 1 + 2/n 

77 


Nil 

Similarity energy 
Dilational energy 



Finally, we remark that conserved densities T with the property 



n ( ) ^ ,4.20) 

are known as a Hamiltonian Casimir. Such conserved densities are distinguished by having no 
correspondence to any symmetry of the Euler equations (I2.ip - (l2.3p . Our classification shows 
that the only vorticity Gasimirs and kinematic Gasimirs admitted by the fluid Hamiltonian 
(14. ip consist of the conserved densities for helicity (13.70 . enstrophy (13.90 . and mass (12.150 . 
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5. Summary and Concluding Remarks 



For the Euler equations fl2.1l) - fl2.3p governing compressible isentropic fluid flow in n > 1 
dimensions, we have directly classified all nontrivial kinematic and vorticity conservation 
laws f|2.4p by solving the determining equations fl2.8p for conserved densities of the respective 
forms fl2.13p and (13.40 . Alternatively, our classification of conservation laws can be carried 
out by means of multipliers [21 El IH [H] ■ This approach is most easily presented as follows 
using index notation (cf. section 2.2). 

Let T(t,x^ , p,u'', pj,u''j, . . .) be a nontrivial conserved density and let X\t,x^ , p,u'', 
Pj,u^ J, . . .) be a spatial flux vector, given by some functions of the time and space co- 
ordinates t and x\ the fluid density p and fluid velocity u\ and their spatial derivatives p,i, 

^i, etc. with respect to x\ which satisfy a local continuity equation 

VtT + DiX' = (5.1) 

where 

Vt = dt- {pu'),dp - {u=u\, + p~'P'{p)p:)d^. + ■■■ (5.2) 

is the time derivative defined through the Euler equations (I2.25p . If we express (15. ID in terms 
of the total derivative Dt = dt + ptdp + u\dui + ■ ■ ■ , then we obtain an equivalent equation 

DtT + DiX' = {pt + ipu'),i)Q + {ui + u^u\, + p~'P\p)p:)Qi (5.3) 

holding for p and given by arbitrary functions of t and x*, where 

Q = Ep{T), Q, = E^^{T) (5.4) 

are functions of t, x^ , p, u^, pj, u''j etc., and where X* differs from X* by terms that are 
linear homogeneous in the Euler equations (I2.25P and total spatial derivatives of (I2.25p . This 
equation (15.30 is called the characteristic form of the conservation law (15. ip . It establishes, 
firstly, that every nontrivial local conservation law of the Euler equations (12.250 arises from 
multipliers (15. 4p . Secondly, since the spatial Euler operators Ep and E^i annihilate diver- 
gences Z)j9* for any vector function 9*(t, x^ ,p,u'' ,pj,u^ j, . . .), the relation (15.40 shows that 
any two local conservation laws differing by a trivial conserved density of the form DjG* 
have the same multipliers. Thus, there is a one-to-one correspondence between nontrivial 
conserved densities (modulo spatial divergences) and non-zero multipliers. 

Necessary and sufficient equations for determining multipliers [H [11] are given by applying 
variational derivative operators S/6p and 6/6u^ to the characteristic equation (15.30 . yielding 
a linear homogeneous polynomial system in pt, u\, ptj, ulj, etc. whose coefficients must 
separately vanish. The resulting determining equations for Q(t,x^ ,p,u^ ,pj,u^ j,. . .) and Qi(t, 



x^ ,p,u^ ,Pj,u'' J,. . .) consist of 

VtQ + u'D.Q + f'{p)D'Q. = 0, (5.5) 

VtQi + uWjQi + 2u{i,Qi] + pD'Q = 0, (5.6) 

^uiQi — ^u'Qj^ ^pQ — l~-pQ-, ^^u^Q — l~-pQ,i-, C,pQi — ^uiQ^ (^•'^) 



where £„i and Cp denote linearization operators (Frechet derivatives) with respect to and 
p; and £* denote the adjoint linearization operators [H [11]. We now note, first, (15.70 
provides the necessary and sufficient conditions [211 [H] for Q and Qi to have the form of 
variational derivatives of some function with respect to p and m*. (Moreover, this function can 
be constructed explicitly from Q and by means of homotopy integral formulas [^ [H [TT| or 
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by an algebraic scaling formula [T] based on invariance of the Euler equations under dilations 
t \t, X* — i> \x\) Second, we note fl5.5p and (15.61) constitute the adjoint of the determining 
equations for symmetries 

Vtf] + D,{u'fi + pf) = 0, (5.8) 
Vtv' + uW.fi' + u\,f]^ + D,{f'{p)fi) = 0, (5.9) 

where 

X = flit, x^, p, M^ pj, u\j, . . .)dp + fi\t, x\p, M^ p,^, m^^-, . . .)9„. (5.10) 

is the symmetry generator in characteristic form. Thus, multipliers can be characterized 
as adjoint-symmetries that have a variational form [U [TT]. In particular, this formulation 
reduces the determination of multipliers and hence of conservation laws to an adjoint version 
of the determination of symmetries. 

We now list the multipliers for the kinematic conservation laws (I2.15p - (l2.23p and vorticity 
conservation laws (I3.7p - (l3.9p in the following two tables. 



For general equations of state (12.30 : 



Conserved density T 


Description 


Q = ST/6p 


Qi = 6T/6u' 


P 

k 

pu 

p{u^x^ — u^x^) 
p{tu^ — x'^) 
\pu^Uk + p J p''^P{p)dp 

U^Wk 

pfi^/p) 


Mass 
Momentum 
Angular momentum 
Galilean momentum 

Energy 
Helicity (n = 2m + 1) 
Enstrophy (n = 2m) 


1 

u'^ 

u^x^ — u^x^ 
tu^ — x^ 
^u'^Uk + / p-'P'{p)dp 


f{zu/p) - p-^wf\w/ p) 




p{x^6j - x^6f) 
pt6f 
pui 
(m + l)'Wi 
mwijf"{w/ p)p^^ 



For the distinguished polytropic equation of state (12.210 : 



Conserved density T 


Description 


Q = ST/6p 


Qi = 6T/6u' 


^ ]_ riK TP 
^PU Uk + —p " — -C/ 

tE — ^pu^Xk 
t^E - p{tu^ - \x^)xk 


Polytropic energy 
Similarity energy 
Dilational energy 


\u^Uk + (1 + f )/«p^ 
\u^{tuk - Xk) + t{l + f )kP" 
\{tuk - Xk){tu^ - x^) + t2(i + D^pi 


pUi 
pitUi - |xj) 
tp{tUi - Xi) 



The adjoint relation between multipliers and symmetries can be expressed in an explicit 
form through the Hamiltonian formulation of the Euler equations (14. ip : 



where 7i is the Hamiltonian operator. As shown in section 4, the mapping defined by (15.110 
annihilates the multipliers for Hamiltonian Casimirs consisting of the conserved densities for 
mass, helicity and enstrophy. From the specific form of the multipliers for the remaining con- 
served densities — momentum, angular momentum, Galilean momentum, energy, similarity 
energy and dilational energy — given in the preceding two tables, we find 

Q = r[\u'ui + y p-^P\p)dp) - iiu' + a, Qi = pirm - 6) (5.12) 

and 

-f]= (ep),i + rpu -ft = e.-w'/' + + rul - a; = {^,u^ - a)/ + 2^juj^' + rul (5.13) 
where ^j, r, a are shown in the following table. 
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T 


a 


Symmetry 


Conserved density 


Equation of State 


-Si 








Space translations 


Momentum 


general 


-25 fx''^ 








Rotations 


Angular momentum 


)) 









Galilean boosts 


Galilean momentum 


n 





1 





Time translation 


Energy 


J) 




t 





Similarity scaling 


Similarity energy 


polytropic 7 = 1 + 2/n 


tXi 


e 


1_ rpj rif . 
2 •A' J 


Galilean dilation scaling 


Dilational energy 





It is straightforward to show that ^j, r, a satisfy the following system of equations: 



2^(i,i) = '^tSij, ^it = <7,i, <Jt = 0, Ti = 0. (5.14) 

In particular, we note 



^i = (i + tXi, r = c+ {2/n)t(C + {l/n)t\iS 




(5.15) 



with c = const., where 

= l^^ij, ^ = const. (5.16) 
is the equation defining dilational Killing vectors Ci{x^) on M", and 

X[i,j] = 0, X{i,j) = l^^ij, ^ = const. (5.17) 

are the equations defining irrotational dilational Killing vectors Xi{^'^) As a result, the 

symmetries corresponding to the multipliers for the non-Casimir conserved densities under 
the Hamiltonian mapping fl5.12p - fl5.13p have the form of geometrical point symmetries 

X = {{2/n) j e„dt)dt + CO.. + e,^pdp + {{l/n)e,^u' + i^^Si - i\)du^ 

given in terms of ^\t,x^) through (15. 150 - 05. 170 . 
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